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Most low-mass stars form in stellar clusters that also contain massive stars, which are sources of far- 
ultraviolet (FUV) radiation. Theoretical models predict that this FUV radiation produces photodissociation 
regions (PDRs) on the surfaces of protoplanetary disks around low-mass stars, which affects planet 
formation within the disks. We report James Webb Space Telescope and Atacama Large Millimeter 

Array observations of a FUV-irradiated protoplanetary disk in the Orion Nebula. Emission lines are detected 
from the PDR; modeling their kinematics and excitation allowed us to constrain the physical conditions 
within the gas. We quantified the mass-loss rate induced by the FUV irradiation and found that it is 
sufficient to remove gas from the disk in less than a million years. This is rapid enough to affect giant 


planet formation in the disk. 


oung low-mass stars are surrounded by 

protoplanetary disks of gas and dust, 

which have lifetimes of a few million 

years (1-3) and are the sites of planet 

formation (4). Planet formation is lim- 
ited by processes that remove mass from the 
disk, such as photoevaporation (5). This occurs 
when the upper layers of protoplanetary disks 
are heated by x-ray or ultraviolet photons. Ra- 
diative heating increases the gas temperature, 
which brings the local sound speed above the 
escape velocity of the disk, causing the gas to 
escape from the system. Those photons could 
be from the central star (6) or from nearby 
massive stars (7). Because most low-mass stars 
form in clusters that also contain massive stars, 
most protoplanetary disks are exposed to ex- 
ternal radiation, and so they are expected to 
experience photoevaporation driven by ultra- 
violet photons during their lifetime (7-11). 
Theoretical models predict that far-ultraviolet 


(FUV) photons, those with energies below the 
Lyman limit (energy E < 13.6 eV), dominate 
the photoevaporation process. The effect affects 
the disk mass, radius, and lifetime (7, 10, 12-18); 
its chemical evolution (19-21); and the growth 
and migration of any planets forming within 
the disk (22). However, these processes have 
not been directly observed. 

Most observational constraints on the mass- 
loss rates associated with photoionization have 
been obtained for objects in the Orion Nebula 
known as proplyds, in which the ionization of 
FUV-driven photoevaporation flows from disks 
produces comet-shaped ionization fronts (23, 24). 
Modeling of the observed ionization fronts of 
proplyds has indicated mass-loss rates M = 10-8 
to 10 © solar masses per year (Mo year *) (25-27). 
However, those observations did not deter- 
mine the physical conditions (radiation field, 
gas temperature, and density) at the locations 
where the photoevaporation flows are launched. 


forms at the disk surface. Most observational 
tracers of PDR physics (spectral lines of Hg, O, 
and C*) are in the near- and far-infrared wave- 
length ranges. The spatial scale of PDRs in ex- 
ternally illuminated disks is a few hundred 
astronomical units (au), which corresponds to 
angular sizes <1 arc sec (") for the closest star- 
forming clusters (72, 29, 30). 


Images of a photoevaporation flow 


Figure 1 shows optical and near-infrared im- 
ages of the Orion Bar, a ridge in the Orion 
Nebula (37) situated about 0.25 pc southeast 
of the Trapezium Cluster of massive stars. 
The western edge of the bar constitutes the 
ionization front (Fig. 1B), which separates 
regions where the gas is fully ionized and at 
temperature T ~ 10* K from the neutral 
atomic region at T ~ 500 to 1000 K. We in- 
vestigated the source [BOM2000] d203-506 
(hereafter d203-506) (32, 33), a protoplanetary 
disk seen in absorption against the bright 
background, which is located at the following 
coordinates: right ascension 5°35™20°.357 and 
declination —5°25'05".81 (J2000 equinox). Pre- 
vious observations of d203-506 found no sign 
of an ionization front (32-34), indicating that 
the radiation field reaching the disk is dom- 
inated by FUV photons. 

We obtained near-infrared and submillimeter 
observations of d203-506, with the James Webb 
Space Telescope (JWST) and the Atacama Large 
Millimeter Array (ALMA), respectively, both 
at angular resolution ~0.1” (corresponding to 
~40 au at the distance of the Orion Nebula). 
JWST images were taken in multiple broad and 
narrow band filters using the Near-Infrared 
Camera (NIRCam) instrument (35). We also ob- 
tained near-infrared spectroscopic observa- 
tions using the integral field unit (IFU) of the 
Near-Infrared Spectrograph (NIRSpec) instru- 
ment on JWST (35). The ALMA interferometric 
data cubes provided maps of rotational emis- 
sion lines from the molecules HCN and HCO", 
with a velocity resolution of 0.2 km s™ (35). 

Figure 2 compares the JWST and ALMA im- 
ages to archival optical images from the Hubble 
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Space Telescope (HST). The nearly edge-on (35) 
dusty disk of d203-506 is visible in absorption 
in all the HST and JWST images (Fig. 2, A to F) 
but in emission in the 344-GHz (870 um) dust 
continuum (Fig. 2G). It is also seen in emission 
with ALMA in Fig. 2H, which shows the HCN 


Fig. 1. Optical and near-infared images of the 
Orion Bar region. (A) HST optical image centered at 
coordinates right ascension 5°35"20°.183 and 
declination -5°25'06".14. [O ui] filter at 502 nm is in 
blue, Ha filter at 656 nm is in green, and [N 11] filter 
at 658 nm is in red. [Credit: NASA/STScI/Rice 
University/C. O'Dell et al. (47)]. (B) JWST near-infrared 
image of the same region at the same scale. Filters 
centered at 1.4 and 2.0 um are in blue; at 2.77, 3.00, 
3.23, 3.35, and 3.32 um in green; at 4.05 um in orange; 
and at 4.44, 4.80, and 4.70 um in red. In (A) and 
(B), the Orion Bar separates regions where the gas is 
fully ionized (upper right) from those where it is in 
neutral form (lower left). (C) Zoomed-in view of the 
d203-506 disk. Red is the NIRCam image in the 
2.12-um filter, which traces molecular hydrogen; blue is 
the 1.64-um filter, which traces [Fe 1] emission lines; 
and green is the emission in the 1.40-um broad-band 
filter, which traces scattered light. 


Fig. 2. Multiwavelength A 


HST 0.65 um Ha 


(v = 0, J = 4 — 3) line, where v and J denote 
the vibrational and rotational quantum num- 
bers, respectively, at 354.505 GHz (845.664 um), 
which traces cold molecular gas. Figure 2, I 
and E, shows emission maps of ALMA HCO* 
(v = 0, J = 4 — 8) at 356.734 GHz (840.381 um) 


images of the d203-506 
disk. (A) HST optical image 
in a Ha filter (23). (B to 

F) JWST near-infrared 
images (35). [(E) is repro- 
duced with permission from 
(33)]. (G to I) ALMA sub- 
millimeter images (35). In all 
panels, the white-filled 
ellipse indicates the size 


and NIRCam H; (v = 1— 0, J = 3 — 1) at 2.12 um, 
respectively, which both show emission from 
the PDR surrounding the disk and absorption 
at the center. Both the H, rovibrational and 
HCO* rotational emission lines trace warm 
(gas kinetic temperatures T;,, ~ 500 to 1000 K) 


C  d203-506 


JWST 1.64 um [Fe II] x10? 


-10 
x10 8.0 


Jy pixel“? 


and shape of the point 
spread function or recon- 


structed telescope beam, 
and the horizontal bar is 
100 au. The white-dashed 
ellipse in (D) indicates the 
shape of the aperture used 
for the extraction of the 
NIRSpec spectrum in Fig. 4. 
The wavelength and 
physical assignment of 
each image is labeled above 
each panel: in (H) and (I), G 
(4-3) is an abbreviation of 
(v=0,J=4— 3). 1Jy= 
1 x 10° Wm? Hz7. 
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molecular gas in PDRs (37). The PDR is also 
bright in the 3.35-um NIRCam filter (Fig. 2F), 
which is dominated by aromatic infrared band 
(AIB) emission from ultraviolet-excited poly- 
cyclic aromatic hydrocarbon (PAH) molecules. 
PAHs are known to be tracers of PDRs (36) 
and have been previously mapped in a proplyd 
in the Orion Nebula (37). The PDR in d203- 
506 is spatially resolved and extends south 
from the disk, forming a lobe shape. A jet is 
visible in the NIRCam [Fe 1] filter at 1.62 um 
(Fig. 2C). A bright emission spot is present in 
the H, and HCO* images in the northwest- 
ern part of the PDR and is also visible in the 
broad-band filter at 1.4 um (Fig. 2B). The loca- 


Fig. 3. Schematic diagram of 
our interpretation of d203-506. 
Dark brown is the edge-on disk of 
cold molecular gas, which appears 
in absorption in the NIRCam 
images but in emission in the 
ALMA HCN and dust continuum 
images (Fig. 2). Brown arrows 
indicate molecular gas escaping 
from this disk, which feeds the 
photoevaporation flow. This 
produces an envelope around the 
disk (light brown shading), which 
is delimited by the dissociation 
front (orange band), where 
molecular hydrogen is dissociated 
into hydrogen atoms by FUV 
photons (pink arrows) coming 
from the Trapezium Cluster. The 
transition from molecular gas in 
the disk to atomic gas under 


Edge-on disk 


~100 au 


tion of this bright spot coincides with the re- 
gion of interaction between the jet and the 
PDR, which is visible only on the side facing 
the Trapezium Cluster. There is also AIB emis- 
sion in the 3.35-um NIRCam filter at this loca- 
tion (Fig. 2F), which indicates ultraviolet 
excitation. Figure 3 shows a schematic dia- 
gram of our interpretation of the morpholog- 
ical features in d203-506. 


Physical properties of the PDR 


Figure 4 shows the NIRSpec spectrum of d203- 
506 (35). Numerous rovibrational emission lines 
of CO (v =1— Oandv=2-—1),OH(w=1—0), 
CH* (v = 1— 0), and H, (up to J = 15) are 


Bright spot 


ee Lobe 


ultraviolet irradiation constitutes the PDR. Blue shows the jet from the central star, which corresponds to the [Fe 1] 
emission (Fig. 2C). The jet interacts with the envelope, producing a bright emission spot (yellow). The 
surroundings of d203-506 (gray) consist of diffuse atomic gas. 
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Fig. 4. JWST NIRSpec spectrum of d203-506. Colored tick marks indicate detected species, as labeled. 
The broad emission bands at 3.3 and 3.4 um are C-H vibrational emission from PAH molecules. Other 
unlabeled lines are mostly atomic (e.g., [O |] or [Fe 1!]). There are no data at wavelengths 2.40 to 2.50 um and 
4.05 to 4.18 um owing to gaps in the NIRSpec detectors. 
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detected. We interpreted these lines as com- 
ing from the PDR, so we traced the physical 
conditions of gas in that region. We modeled 
(35) the H; lines using the Meudon PDR code 
(38), which calculates the H, excitation in 
PDRs (Fig. 5). We derived the hydrogen num- 
ber density ny and temperature of the gas in 
the H, emitting layer (compare with Figure 
2E) as ny = 5.5 x 10° to 1.0 x 10’ cm™ and 
Teas = 1240 to 1260 K. 

We fitted a Keplerian orbit model to the 
HCN emission map (fig. S2) and used it to 
set an upper limit on the mass of d203-506’s 
central star M, < 0.3 Mo (35). Taking Tyas ~ 
1250 K as determined above, the speed of 
sound ¢s= s = 3.3 km s~, where kz is the 
Boltzmann constant, mm is the mass of the 
hydrogen atom, and ų is the ratio of total mass 
to hydrogen mass of the gas, which we assume 
equals the interstellar value [u = 1.4 (39)]. 

This value of cs exceeds the escape velocity 
at distances from the central star above a crit- 
ical value, defined as the gravitational radius 
t= Ty (40), where Gis Newton’s gravitational 
constant. For M, < 0.3 Mo and Toas ~ 1250 K, 
we found r, < 26 au. This is much smaller than 
the observed radial extent of the H, emission 
Ty, = 132413 au [and its height hy, = 56+13 au 
(35)]. Therefore, the gas in this layer is not 
gravitationally bound and flows outward from 
the disk, at roughly the speed of sound. The 
associated mass flux through the PDR is 7 = 
UMgnycs, and the total mass loss rate isM = 
j x S, where S is the surface area of the H3- 
emitting layer (35). Including the uncertainties 
on TH, Au,» ne, and Tyas (table S1), we calcu- 
lated M = 14 x 10” to 4.6 x 10° Mz year’. 

We also investigated an alternative method 
of determining the mass-loss rate, using 1D 
dynamical models, and found a consistent 
value of M (35). All the physical quantities 
that we derived are listed in table S1. 


Implications for planet formation 


Gas in protoplanetary disks is the raw material 
from which giant planets form, so mass loss due 
to photoevaporation can limit the formation of 
such planets. The effects of FUV radiation depend 
on the stellar mass, which sets the strength of the 
gravitational field that acts to retain the gas. 
Theoretical models of planet growth under 
the influence of external FUV photoevapora- 
tion have predicted that FUV radiation fields 
with intensity Go = 500 [where Go is the ratio of 
the ultraviolet radiation field to the average 
value in the interstellar medium (47)] sup- 
press giant planet formation around stars 
with masses < 0.5 Mọ (22). Our results for d203- 
506 are consistent with this prediction: We 
found that M, < 0.3 Mo and the radiation 
field is Go $10° (35), whereas the mass-loss 
rate M = 1.4 x 107’ to 4.6 x 10°° Mo year ' 
implies a disk depletion timescale t= Ma/M < 
0.13 million years, where Mg is the disk mass 
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Fig. 5. Observed and modeled Hp. Line intensities for d203-506. Blue squares show the observed line 
intensities (table S2), with error bars enclosing 50% uncertainty, estimated from the x? (35). The 
instrumental uncertainties are smaller than the symbols. Colored lines are the best-fitting models (35) 
using small dust grains (orange) or large dust grains (green); circles show the model intensities of each 
line. H2 line identifications are abbreviated; the full quantum levels corresponding to this notation are 


listed in table S2. 


(35). This is faster than even very early planet 
formation (42, 43). 

A positive correlation has been observed 
between stellar mass and the frequency of 
Jupiter-mass exoplanets orbiting those stars 
(44, 45), which we suggest could be due to 
FUV radiation in stellar clusters during the 
planet formation process. Dynamical and com- 
positional studies of the Solar System indicate 
that it formed in a stellar cluster that contained 
one or more massive stars (46), and so it might 
have been affected by FUV radiation. 
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